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L-SAOMs  for  Relational  Events 
^ — SAOMs  &  Relational  Events 


I —  Introduction 

A  convenient  and  powerful  combination... 


h  Using  relational  event  datasets  to  study  social  network 
dynamics  can  have  many  practical  up-sides. 

a  SAOM  effects  (which  are  based  on  digraphs)  are  convenient 
(evidenced  by  the  large  amount  of  effects  that  have  already 
been  studied). 
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^ — SAOMs  &  Relational  Events 


I — Practical  considerations 

Practical  considerations  for  relational  event  data 


h  Collecting  high  quality  relational  event  datasets  can  be  cheap 
and  easy. 
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I — Practical  considerations 

Practical  considerations  for  relational  event  data 


a  Collecting  high  quality  relational  event  datasets  can  be  cheap 
and  easy. 

a  Datasets  can  be  quite  large  (and  correspondingly  contain  lots 
of  information). 
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I — Practical  considerations 

Practical  considerations  for  relational  event  data 


a  Collecting  high  quality  relational  event  datasets  can  be  cheap 
and  easy. 

a  Datasets  can  be  quite  large  (and  correspondingly  contain  lots 
of  information). 

a  We  can  “follow  the  sound  of  marching  feet.” 
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L-SAOMs  for  Relational  Events 
^ — SAOMs  &  Relational  Events 


I — Practical  considerations 

Practical  considerations  for  “relational  event  history 
models” 


Butts  [2008]  proposes  to  study  relational  event  histories  with  a  Cox 
model. 
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^ — SAOMs  &  Relational  Events 


I — Practical  considerations 

Practical  considerations  for  “relational  event  history 
models" 


Butts  [2008]  proposes  to  study  relational  event  histories  with  a  Cox 
model. 

Devising  good  effects  with  relational  event  histories  can  be  less 
convenient  than  with  SAOMs: 

■  How  do  we  handle  the  passing  of  time?  How  much  more 
important  is  something  that  happened  last  week  than  last 
year? 
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I — Practical  considerations 

Practical  considerations  for  “relational  event  history 
models" 


Butts  [2008]  proposes  to  study  relational  event  histories  with  a  Cox 
model. 

Devising  good  effects  with  relational  event  histories  can  be  less 
convenient  than  with  SAOMs: 

■  How  do  we  handle  the  passing  of  time?  How  much  more 
important  is  something  that  happened  last  week  than  last 
year? 

a  How  do  we  model  “second  order  terms”  like  transitivity, 
balance,  and  betweenness? 


<  □ 


3 


>0  0.0 


L-SAOMs  for  Relational  Events 
^ — SAOMs  &  Relational  Events 


Idea:  model  the  evolution  of  the  affective  network  with 

SAOMs 


...but  where  are  the  digraphs? 

n  A  fairly  common  practice  is  to  do  some  "binning”  to  make 
data  that  looks  like  digraph  panel  data 

■  Intuitively,  we  suppose  that  more  interaction  in  a  time  interval 
signals  an  affective  tie. 


□ 


3 


>0  0.0 


L-SAOMs  for  Relational  Events 
^ — SAOMs  &i  Relational  Events 
I —  Binning 


We  can  easily  construct  weighted  networks  over  some  time 
intervals... 
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L-SAOMs  for  Relational  Events 
^ — SAOMs  &i  Relational  Events 
I —  Binning 


And  use  some  "threshold”  to  yield  a  digraph: 
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^ — SAOMs  &i  Relational  Events 


Here's  a  higher  threshold: 
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^ — SAOMs  &  Relational  Events 


I —  Binning 

But  really  we’ve  traded  one  set  of  difficulties  for  another: 


■  How  much  interaction  constitutes  a  tie? 
a  How  do  you  choose  intervals? 

a  How  do  we  minimize  information  loss? 

■  How  can  we  be  sure  that  artifacts  don’t  creep  into  our  results? 
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L-SAOMs  for  Relational  Events 
I —  Data 

^ —  Background 

Data:  The  Ikenet  Study 


a  Kate  Coronges,  in  collaboration  with  the  USMA  Network 
Science  Center,  collected  the  email  traffic  among  22 
mid-career  Army  officers  over  a  13-month  period 

b  We  will  use  one  covariate  which  indicates  whether  the  officer 
is  a  captain  (lesser  rank)  or  a  major  (greater  rank). 

■  For  emails  that  are  not  (a)  broadcast  emails  (i.e.  those  emails 
sent  to  the  whole  group)  or  (b)  self-sent  emails,  we  add  one 
“relational  event”  per  recipient  in  a  random  order.  This 
process  yields  8819  events. 
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I —  Data 

^ —  Background 

Data:  The  Ikenet  Study 


a  Kate  Coronges,  in  collaboration  with  the  USMA  Network 
Science  Center,  collected  the  email  traffic  among  22 
mid-career  Army  officers  over  a  13-month  period 

■  We  will  use  one  covariate  which  indicates  whether  the  officer 
is  a  captain  (lesser  rank)  or  a  major  (greater  rank). 

■  For  emails  that  are  not  (a)  broadcast  emails  (i.e.  those  emails 
sent  to  the  whole  group)  or  (b)  self-sent  emails,  we  add  one 
“relational  event”  per  recipient  in  a  random  order.  This 
process  yields  8819  events. 

a  Bonus:  friendship  surveys  are  collected  at  each  month 
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I —  Descriptives 


Frequency  of  Emails  by  Date 
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I —  Descriptives 


Email  Frequency  by  Hour 
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I —  Descriptives 
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I — Ikenet  Analysis  with  Binning 
I — Some  a  priori  decisions 

Binning 


It's  not  obvious  how  we  determi 
■  time  intervals 
0  dichotimization  thresholds 
to  create  waves  out  of  this  data 
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I — Some  a  priori  decisions 

Binning 


It's  not  obvious  how  we  determine 

■  time  intervals 

■  dichotimization  thresholds 
to  create  waves  out  of  this  data. 

To  simplify  the  example,  we  (quite  arbitrarily)  decide  to  create 
intervals  by  month  of  the  year,  yielding  13  digraphs  in  the  panel. 
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I — Some  a  priori  decisions 

Binning 


It's  not  obvious  how  we  determine 

■  time  intervals 

■  dichotimization  thresholds 
to  create  waves  out  of  this  data. 

To  simplify  the  example,  we  (quite  arbitrarily)  decide  to  create 
intervals  by  month  of  the  year,  yielding  13  digraphs  in  the  panel. 
How  do  various  threshold  values  affect  these  digraphs? 
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I — Ikenet  Analysis  with  Binning 
I — Thresholds 
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I — Ikenet  Analysis  with  Binning 
L-  Model 

Model 


Since  there  is  no  clear  choice  of  a  threshold  at  which  the  network 
density  is  uniquely  plausible,  we’ll  run  estimations  of  the  following 
model  at  various  thresholds: 
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Model 


Since  there  is  no  clear  choice  of  a  threshold  at  which  the  network 
density  is  uniquely  plausible,  we’ll  run  estimations  of  the  following 
model  at  various  thresholds: 

■  Outdegree:  “do  I  like  to  have  ties?” 

hi  Reciprocity:  “do  I  like  to  have  ties  to  match  incoming  ties?” 
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Model 


Since  there  is  no  clear  choice  of  a  threshold  at  which  the  network 
density  is  uniquely  plausible,  we’ll  run  estimations  of  the  following 
model  at  various  thresholds: 

■  Outdegree:  “do  I  like  to  have  ties?” 

a  Reciprocity:  “do  I  like  to  have  ties  to  match  incoming  ties?” 

■  Trans.  Triplets:  “do  I  like  to  have  ties  to  ties-of-ties  ?” 

■  Balance:  “do  I  like  my  ties  to  match  my  current  alters’  ties?” 
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Model 


Since  there  is  no  clear  choice  of  a  threshold  at  which  the  network 
density  is  uniquely  plausible,  we’ll  run  estimations  of  the  following 
model  at  various  thresholds: 

■  Outdegree:  “do  I  like  to  have  ties?” 

a  Reciprocity:  “do  I  like  to  have  ties  to  match  incoming  ties?” 

■  Trans.  Triplets:  “do  I  like  to  have  ties  to  ties-of-ties  ?” 

■  Balance:  “do  I  like  my  ties  to  match  my  current  alters’  ties?” 

s  In.  Popularity:  “do  I  like  to  tie  to  alters  with  lots  of  incoming 
ties?” 

a  In.  Activity:  “do  I  like  to  create  more  ties  when  I  recieve  lots 
of  ties?" 
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Model 


Since  there  is  no  clear  choice  of  a  threshold  at  which  the  network 
density  is  uniquely  plausible,  we’ll  run  estimations  of  the  following 
model  at  various  thresholds: 

■  Outdegree:  “do  I  like  to  have  ties?” 

a  Reciprocity:  “do  I  like  to  have  ties  to  match  incoming  ties?” 

■  Trans.  Triplets:  “do  I  like  to  have  ties  to  ties-of-ties  ?” 

a  Balance:  “do  I  like  my  ties  to  match  my  current  alters’  ties?” 

a  In.  Popularity:  “do  I  like  to  tie  to  alters  with  lots  of  incoming 
ties?” 

a  In.  Activity:  “do  I  like  to  create  more  ties  when  I  recieve  lots 
of  ties?" 

■  Same  Rank:  “do  I  like  to  create  ties  with  other  officers  of  the 
same  rank?” 
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I — Ikenet  Analysis  with  Binning 
I— Results 

MoM  Estimation  [Snijders,  2001]  Results 
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0.11 
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Balance 
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In.  Popularity 

0.13 
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L-SAOMs  for  Relational  Events 
I — Ikenet  Analysis  with  Binning 
I — Summary 

Summary  on  binning 


We  get  convenient  SAOM-based  inference  (easy  to  interpret 
effects)  out  of  relational  event  data,  but... 

a  Choice  of  binning  parameters  is  not  obvious  and  may  have 
important  consequences  on  results 

■  How  much  information  was  lost  by  binning? 

■  Did  we  introduce  artifacts? 
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I — L-SAOMs  for  Relational  Events 
I —  Introduction 

Bringing  relational  events  into  the  SAOM 


Let's  consider  how  we  might  bring  the  relational  events  into  the 
SAOM: 

a  We  consider  the  relational  event  mode  as  just  another  aspect 
over  which  ego  has  complete  control. 
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Let's  consider  how  we  might  bring  the  relational  events  into  the 
SAOM: 

a  We  consider  the  relational  event  mode  as  just  another  aspect 
over  which  ego  has  complete  control. 

■  An  ego  may  then  choose  to  modify  her  outgoing  ties,  or  to 
send  a  relational  event  to  an  alter. 
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Bringing  relational  events  into  the  SAOM 


Let's  consider  how  we  might  bring  the  relational  events  into  the 
SAOM: 

a  We  consider  the  relational  event  mode  as  just  another  aspect 
over  which  ego  has  complete  control. 

s  An  ego  may  then  choose  to  modify  her  outgoing  ties,  or  to 
send  a  relational  event  to  an  alter. 

■  As  the  network  is  unobserved,  this  is  a  hidden  Markov  model 
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I — L-SAOMs  for  Relational  Events 
I — Notation 


We  can  think  of  this  process  in  terms  of  ordered  ministep  tuples 
v  =  (i,j,  k,t)eV  where: 
m  i  e  Af  indexes  the  ego 
o  j  E  J\T  indexes  the  alter 
a  k  G  /C  indexes  the  aspect 
■  t  E  indexes  time 

where  Af  is  the  actor  set,  /C  is  the  set  of  networks  and  modes. 
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We  can  think  of  this  process  in  terms  of  ordered  ministep  tuples 
v  =  M)fV  where: 
m  i  e  Af  indexes  the  ego 
m  j  e  Af  indexes  the  alter 
a  k  G  /C  indexes  the  aspect 
■  t  G  SK  indexes  time 

where  Af  is  the  actor  set,  K,  is  the  set  of  networks  and  modes. 

For  convenience,  these  indices  can  be  used  as  functions,  e.g.  i(v 5) 
refers  to  the  ego  of  the  5-th  ministep. 
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We  can  think  of  this  process  in  terms  of  ordered  ministep  tuples 
v  =  M)fV  where: 
m  i  e  Af  indexes  the  ego 
m  j  e  Af  indexes  the  alter 
a  k  G  /C  indexes  the  aspect 
■  t  G  SK  indexes  time 

where  Af  is  the  actor  set,  K,  is  the  set  of  networks  and  modes. 

For  convenience,  these  indices  can  be  used  as  functions,  e.g.  i(v 5) 
refers  to  the  ego  of  the  5-th  ministep. 

Also  for  convenience,  let  Vr  =  {v^  G  V  :  t(Vfa)  <  r} 
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I — L-SAOMs  for  Relational  Events 
L-  Model 


If  we  let 

■  po(X*)  =  po(X*\6)  represent  the  PMF  of  the  initial  network 

0  fv{vb )  =  f\z(vb|Vt(Vfa),x*,  6)  represent  the  conditional  PDF 
given  all  ministeps  occurring  before  v 5  and  the  initial  network. 
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If  we  let 

■  po(X*)  =  po(X*\9)  represent  the  PMF  of  the  initial  network 

0  fv{vb)  =  f\z(vb|Vt(Vfa),x*,  6)  represent  the  conditional  PDF 
given  all  ministeps  occurring  before  v b  and  the  initial  network. 

then  the  likelihood  is 

L{6\V ,  x*)  =  p0{x*)  fv(vb)  (1) 

vb£V 
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For  po(x*),  we  are  free  to  choose  anything.  In  the  context  of  the 
Ikenet  data,  we  use  a  simple,  dyad-independent  model: 

Po(X)=  1]  p(Xij  =  xiJ,Xji  =  Xji\0)  (2) 
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For  po(x*),  we  are  free  to  choose  anything.  In  the  context  of  the 
Ikenet  data,  we  use  a  simple,  dyad-independent  model: 

Po(X)=  1]  p(Xij  =  xiJ,Xji  =  Xji\0)  (2) 


We  include  terms  in  p(X,y  =  1| xy;,$)  for  density  and  reciprocity. 
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For  fv(vb),  we  decompose  the  joing  wait-time/ego/aspect  density, 
and  alter  selection  density: 


fv(yb)  =  fikt{i{vb),  k(vb),  t{vb)}  fj\ikt{j{vb)} 
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For  fv(vb),  we  decompose  the  joing  wait-time/ego/aspect  density, 
and  alter  selection  density: 


fv(yb)  =  fikt{i{vb),  k(vb),  t(vb)}  fj\ikt{j{vb )} 


Just  as  in  the  regular  SAOM,  we  suppose  f;kt  is  a  negative 
exponential  process,  and  we  have  a  different  rate  for  each  aspect 
K. 
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For  a  network  aspect  k,  our  alter  selection  density  is  the  SAOM 
alter  selection  density: 


g{i  ~*j) 

En&Mg{l  ^  k) 


<  &  ►  «  1  ►  «  1  ► 
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For  a  relational  event  mode  aspect  k,  our  alter  selection  density 
can  be  anything,  but  we  choose  a  parsimonious  “observation 
function” 


r_i  _  exp{7Xy(t)} 
miJ 1  En^wexp{7x/n(t)}  ' 


■*  □  ► 


<  &  ► 


'OQ ,o 


L-SAOMs  for  Relational  Events 
I — L-SAOMs  for  Relational  Events 
L-  Model 


For  a  relational  event  mode  aspect  k,  our  alter  selection  density 
can  be  anything,  but  we  choose  a  parsimonious  “observation 
function” 


r_i  _  exp{7Xy(t)} 
miJ 1  En^wexp{7x/n(t)}  ' 


The  “network  effect”  7  represents  the  tendency  for  an  ego  to  want 
to  send  relational  events  to  alters  with  whom  they  have  a  network 
tie. 
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DAG  for  the  L-SAOM 


x(r(ei)) - *-  x{Li) 


Time  t  =  r(ei) 
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In  principle,  estimation  is  a  straightforward  MCMC-MLE  scheme 
[Snijders  et  al. ,  2010]. 
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In  principle,  estimation  is  a  straightforward  MCMC-MLE  scheme 
[Snijders  et  al. ,  2010]. 

We’ve  added  a  few  proposals  to  the  MCMC  scheme  to  account  for 
the  lack  of  so-called  "parity”  conditions.  The  details  are  pretty 
involved. 
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In  principle,  estimation  is  a  straightforward  MCMC-MLE  scheme 
[Snijders  et  al. ,  2010]. 

We’ve  added  a  few  proposals  to  the  MCMC  scheme  to  account  for 
the  lack  of  so-called  "parity”  conditions.  The  details  are  pretty 
involved. 

For  the  sake  of  time,  we’ll  surf  right  over  these  details.  They  are 
implemented  in  a  .NET  port  of  RSiena: 

github . com/ JLospinoso/sie 
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6 

s.e.  9 

6 

s.e.  6 

Network  Dynamics 

Observation 

Outdegree 

-3.19 

0.70 

Network 

2.30 

0.71 

Reciprocity 

0.84 

0.15 

Initial  Network 

Trans.  Triplets 

-0.03 

0.06 

Outdegree 

-3.49 

0.53 

Balance 

0.04 

0.07 

Reciprocity 

7.24 

0.94 

In.  Popularity 

0.07 

0.01 

Pacing 

In.  Activity 

0.07 

0.08 

Network 

4.73 

1.05 

Same  Rank 

0.15 

0.06 

Email 

5.99 

0.27 
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How  do  these  compare  with  the  actual  friendship  survey  data? 
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How  do  these  compare  with  the  actual  friendship  survey  data? 


Email 

Friendship 

6 

s.e.  9 

0 

s.e.  8 

Outdegree 

-3.19 

0.70 

-0.95 

0.20 

Reciprocity 

0.84 

0.15 

1.12 

0.19 

Trans.  Triplets 

-0.03 

0.06 

-0.24 

0.23 

Balance 

0.04 

0.07 

0.21 

0.13 

In.  Popularity 

0.07 

0.01 

0.09 

0.03 

In.  Activity 

0.07 

0.08 

0.03 

0.10 

Same  Rank 

0.15 

0.06 

0.24 

0.09 
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Binning  results  again: 
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Summary 


s  Using  SAOM  effects  to  model  relational  event  data  is  a 
convenient  and  powerful  combination. 

■  Binning  might  not  be  such  a  great  idea. 

■  We  can  extend  SAOMs  to  include  a  relational  event  mode  and 
use  an  observation  function. 

■  In  our  study,  inferential  results  are  very  similar  for  friendship 
and  email  data. 
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